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Abstract 



We describe a novel phenomenon in which vortices are produced due to 
resonant oscillations of a scalar field which is driven by a periodically varying 
temperature T, with T remaining much below the critical temperature Tc- 
Also, in a rapid heating of a localized region to a temperature below Tc, far 
separated vortex and antivortex can form. We compare our results with recent 
models of defect production during reheating after inflation. We also discuss 
possible experimental tests of our predictions of topological defect production 
without ever going through a phase transition. 

PACS numbers: 98.80.Cq, 11.27.+d, 67.40.Vs 
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Recently, a lot of interest has been focussed on the formation and consequences of topo- 
logical defects both in the context of particle physics models of the early universe [0] as 
well as in condensed matter systems [0. In all the investigations, defects are produced 
by two different mechanisms. They could be produced by external influence, such as flux 
tubes in superconductors in an external magnetic field or vortices in superfiuid helium in 
a rotating vessel. The only other method of producing topological defects is either during 
a phase transition (via Kibble mechanism [^), or due to thermal fluctuations (with defect 
density suppressed by the Boltzmann factor). It has been recently demonstrated by some 
of us IQ that defects can be produced via a new mechanism due to strong oscillations, and 
subsequent flipping, of the order parameter (OP) field during a phase transition. 

In this letter we discuss a new phenomenon where defects can form without ever going 
through a phase transition. In our model, defects form due to resonant oscillations of the 
OP field at a temperature which is periodically varying, but remains much below the critical 
temperature Tc. We also show that in rapid heating of a localized region to a temperature 
below Tc, far separated vortex and antivortex (or a large string loop in 3+1 dimensions) can 
form. Our model Lagrangian describes a system with a spontaneously broken global U(l) 
theory in 2+1 dimensions. The Lagrangian is expressed in terms of scaled, dimensionless 
variables. 



Here $ = $i + i$2 is a complex scalar field with magnitude 0. T is the temperature of 
the system (with Tc = 1) and a is a dimensionless parameter. We take a = 4. We get similar 
results with linear temperature dependence in the effective potential. We emphasize that 
the basic physics of our model resides in the time dependence of the effective potential. We 
achieve this by using a time dependent temperature. One could also do this by periodically 
varying some other parameter such as pressure (which may be experimentally more feasible), 
or even a time-dependent external electric or magnetic field (say, for liquid crystals 0). The 
oscillatory temperature dependence in our model leads to field equations which are similar 
to those with an oscillating inflaton field coupled to another scalar field in the models of 
post-inflationary reheating [§-[10| (more precisely, to the case with spontaneously broken 
symmetry for the scalar field [^]), with the oscillating temperature playing the role of the 
inflaton field. Analysis of growth of fluctuations shows the existence of exponentially growing 
modes in this case implying that fluctuations will grow rapidly for modes with wave vector 
lying in the resonance band. However, there are also crucial differences between our results 
and those in the literature [|6|-p!0[l, as we will explain later. 

We use the following equations for field evolution. 



Here V{(j)) is the effective potential in Eqn.(l) and r/ is a dissipation coefficient. We solve 
Eqn. (2) using a second order staggered leapfrog algorithm. We do not include a noise term 
in Eqn. (2). The basic physics we discuss does not depend on noise. Also, as discussed above, 
time dependence of the effective potential could be thought to arise from some other source 
than the temperature (with temperature kept low to suppress any thermal fluctuations). 



c=-id,^^w<^)--i<i^'-ir--T'<i>'. 



(1) 



+ r]d^i/dt - v'$^ + ^'(0) = 0. 



(2) 
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Here we mention that the effect of noise on the growth of fluctuations has been studied 



in ref. in the context of reheating after inflation. It is shown in ref. [|lT| that small 
amplitude noise does not affect growth of fluctuations for unstable modes. In a future work 
we will study the effect of noise in our model, and also explore connections with the well 
studied phenomenon of stochastic resonance in condensed matter systems |T2[ . 

We first discuss a situation in which a central patch in the system is instantaneously 
heated from T = to T = Tq < Tc and maintained at that temperature, while the region 
outside the patch is kept at T = 0. Field is evolved by Eqn.(2), with T being the local 
temperature. The temperature profile of the region separating the central patch from the 
surrounding region is taken to be T(r) = (To/2)(l — tanh{{r — i?)/A)); R is the radius of 
the (circular) central patch and A is the wall thickness. (Our results do not depend on the 
profile used for this boundary region.) We use a 500 x500 lattice with physical size equal 
to 80 X 80. Time step At was taken to be 0.008. We take = 24 and A = 0.32. Initial 
magnitude of $ was taken to be equal to the T = vacuum expectation value (vev), while 
its phase varied linearly from to 27r along the Y-axis, being uniform along the X-axis. 
(This choice facilitated the use of periodic boundary conditions to evolve the field. Defect 
production happens with much smaller phase variation also.) The instantaneous heating of 
the central region with Tg = 0.72 (with rj = 0.4), led to destabilization of $ in that region. 
The field, while trying to relax to its new equilibrium value, overshot the central barrier and 
ended up on the opposite side of the vacuum manifold. As discussed in [Q, this results in 
the formation of a vortex at one boundary of the flipped region, with an anti- vortex forming 
at the opposite boundary, as shown in Fig.l. This occurs because overshooting of $ (in a 
localized region) leads to a discontinuous change in its phase by vr {flipping of $). Note that 
these defects could not have formed via the Kibble mechanism as the system never goes 
through a symmetry breaking phase transition in our case. Similarly, thermal fluctuations 
can not directly produce a pair in which defect and antidefect are so far separated. 

This shows that above scenario can be used to provide an unambiguous experimental 
verification of the flipping mechanism described in for defect production. For example, 
one can locally heat a central portion of superfluid ^He system to a temperature Tq < {Tq 
should be large enough to allow for overshooting ||13[). In 3 + 1 dimension, this would lead to 
the formation of a string loop, with a size of order of the size of the heated central region (in 
2 + 1 dimensions one will get far separated vortex-antivortex pair). The requirement of small 
phase variation can be easily achieved for superfluid helium by allowing for small, uniform, 
superflow across the system. (For superconductors, one can allow small supercurrent.) The 
fact that the formation of such a large string loop does not require a phase transition, makes 
this a clean signal, unpolluted by the presence of smaller string loops. 

We now describe the resonant production of vortex-antivortex pairs, induced by peri- 
odically varying temperature T. T (in Eqn.(l)) is taken to be spatially uniform, with its 
periodic variation given by T{t) = Tq + TaSm{ut). The choice of frequency u was guided 
by the range of frequency required to induce resonance for the case of spatially uniform 
field (evolved by Eqn.(2)). We find that resonance happens when u lies in a certain range. 
(We are assuming that for the relevant range of u here, the system can be considered in 
quasi-equilibrium so that the use of temperature dependent effective potential makes sense.) 
This frequency range for which resonant defect production occurs depends on the average 
temperature {Tq) as well as on the amplitude of oscillation T^, with the range becoming 
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larger as To approaches Tc. 

Initially the magnitude of $ is taken to be equal to the vev at T = Tq, over the entire 
lattice. The phase ^ of $ was chosen to have small random variations at each lattice point, 
with 9 lying between 20° and 40°. We have carried out simulations for other widely different 
initial configurations as well. In one case we took the domain structure developing at late 
stages in above simulations, and used it as the initial data for evolution with a different set 
of parameters. In another case, we took the lattice to consist of only four domains with 
small variations of 6 between domains. Defects are always produced as long as there is some 
non-uniformity in the initial phase distribution. Although we find that at early times defect 
formation depends on the initial configuration, the asymptotic average defect density seems 
to remain unchanged. Note that the analysis based on the Mathieu equation, ||6HlO| in 
the context of inflationary models, suggests that the growth of fluctuations may depend on 
the initial configuration. However, such analysis cannot be trusted for times beyond which 
fluctuations become too large. 

Fig. 2 shows a portion of the lattice for uo = 1.19 and t] = 0.013, with four vortex- 
antivortex pairs. We count only those pairs which are separated by a distance larger than 
2m~^. We emphasize again, formation of these vortices cannot be accounted for by Kibble 
mechanism, or by the thermal production of defects. To compare our results with the defect 
density expected in the case of thermally produced defects, we have estimated the energy of 



a defect-antidefect pair, using the numerical techniques in ref. |jTj]. We find it to be about 
2.5 (for To = 0.31) when the separation between the vortex and the antivortex is 2nijl (with 
temperature dependent mu)- If we take thermally produced defect density to be of order 
~ T^ exp(- impair /T), then even with T = To + Tq = 0.435, the defect density is only about 
6 X 10~^, which is less than the defect densities we find. 

We now discuss differences between defect production in our model and those discussed 



in refs. IpHlOf in the context of reheating after inflation. In these works, non-thermal fluctua- 
tions grow large enough to restore the symmetry Topological defects are produced when 
the symmetry subsequently breaks due to rescattering effects and/or universe expansion. 
However, in our case there is no symmetry restoration ever (see, also, ref. |jl5|). To check 
this we plot the probability distribution of $ in Fig.3. This clearly shows that the peak of 
the probability distribution lies around a valley of radius ^ 0.9 (for To = 0.375) and not at 
the zero of $ as would have been the case if symmetry had been restored 0. In all the cases 
we have studied, we never find the peak of the probability density of $ to be centered near 
$ = (except for early transient stages when $ in the entire lattice flips through $ = 0). 
Further, detailed field plots show that the defect production in our model happens via the 
flipping mechanism and not by the formation of domain structure (via Kibble mechanism) 
as would be expected if defects were produced due to non-thermal symmetry restoration 
and its subsequent breaking. Hence the defects are always produced in pairs in our case. 
Also, the pair production continues for all times because of continued localized flipping of $ 
under the influence of the periodically oscillating temperature (equivalently, an oscillating 
inflaton field). 

We have carried out simulations for various values of To, (with uj = 1.19, T^ = 0.125, 
and f] = 0.013 kept fixed), starting with similar initial configurations. As a function of 
time, defect density rises from zero, and eventually fluctuates about an average value. We 
carry out the evolution till t ~ 1000 when this average defect density becomes reasonably 
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constant. For this set of parameters, the smallest value of Tq, for which we observe resonant 
defect production, was equal to 0.31. (For linear temperature dependence in V{(j)) in Eqn.(l) 
we get defect production at much lower temperatures, with smallest value of Tq = 0.10 and 
Ta = 0.08, uj = 1.19.) Fig.4 shows temporal variation of defect density at large times. The 
fluctuations in defect density are due to pair annihilations and due to pair creation of vortices 
which keeps happening periodically because of localized flipping of $. 

Fig. 5 shows plot of the asymptotic average defect density (for 900 < t < 1000) as a 
function of Tq. One can see that defect density peaks at a value of To ~ 0.5. This drop 
in average defect density for larger temperatures can be due to decrease in rriH. Due to 
flatter Veff near the true vacuum, $ does not gain too much potential energy when Ve// 
changes due to change in temperature (with Ta = 1.0), which may make flipping of $ 
difficult (even though the barrier height is smaller now). A near perfect Gaussian fit to the 
values of defect density in Fig. 5 may be indicative of a possibility of establishing some sort 
of correspondence with a thermodynamic phase of the system with some different effective 
temperature. Investigation of this possibility requires accurate determination of asymptotic 
defect density and its dependence on different parameters, such as Tq, Ta, and u. (As one 
can see from Fig.4, average defect densities show slow variation over a large time scale for 
certain cases.) 

We find that there is an optimum value of frequency uj for which the defect density is 
maximum. It is reasonable to expect that if uj is too large, it would be difficult for the field 
to overshoot the central barrier since larger frequency will lead to averaging out of the force 
on $ due to the rapidly changing effective potential. On the other hand, if u is too small, 
the field may have enough time to relax to its equilibrium value without being sufficiently 
destabilized, making the overshooting of $ more difficult. By keeping all other parameters 
fixed, (Tq = 0.5, Ta = 0.13 and r] = 0.013) we have explored the behavior of defect density 
for different values of frequency. (For the above parameters, resonant defect production 
happens for 0.94 < u; < 1.31). For ou = 0.94, 1.0, 1.13, 1.19, and 1.31, we find the average 
defect density to be 0.0048, 0.0061, 0.0059, 0.0058, and 0.0042 respectively. 

We have made some check on the dependence of the asymptotic average defect density 
on Ta- For Tq = 0.5, u = 1.19 and r] = 0.013, the average defect density is larger for 
Ta = 0.13 compared to the case when Ta = 0.25. (Again, note that we are focusing on 
asymptotic defect density.) When we increase t], it results in suppression of the defect 
density since damping makes it difficult for $ to overshoot the barrier. Moreover, as t] is 
gradually increased, the minimum value of Tq required to induce resonant oscillations also 
increases and for i] > 0.63, resonant production of vortices is not possible for Tq < T^. We 
also find that the range of frequency for which resonant defect production occurs, becomes 
narrower on increasing r]. 

In conclusion, we have demonstrated an interesting phenomenon where topological de- 
fects can form at very low temperatures, without the system ever going through the phase 
transition. It should be easy to experimentally verify this scenario of defect production 
in many condensed matter systems, such as superfiuid helium, or superconductors. There 
are also important implications of our results for defect production during reheating after 
infiation. Earlier studies p|-p!0| have focussed on defect production due to non-thermal sym- 
metry restoration and its subsequent breaking. Our results show the possibility that defects 
may be produced due to oscillations of the infiaton field via the fiipping mechanism even 
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when there is no symmetry restoration ever. This may lead to defect production under more 
generic conditions than discussed in the hterature. Our results suggest the very interesting 
possibility that various models of re-heating after inflation can have analogs in the condensed 
matter systems which may provide possible experimental tests for the predictions of these 
models. 



6 



REFERENCES 



[1] A.Vilenkin and E.P.S.Shellard, "Cosmic Strings and Other Topological Defects", (Cam- 
bridge University Press, Cambridge, 1994). 
[2] W.H. Zurek, Phys. Rep. 276, 177 (1996). 
[3] T.W.B. Kibble, J. Phys. A9, 1387 (1976). 

[4] S. Digal, S. Sengupta, and A.M. Srivastava, Phys. Rev. D55, 3824 (1997); ibzd Phys. 

Rev. D56, 2035(1997). 
[5] S. Digal, S. Sengupta, and A.M. Srivastava, Phys. Rev. D58, 103510 (1998). 
[6] J.H. Traschen and R.H. Brandenberger, Phys. Rev. D42, 2491 (1990); Y. Shtanov, J. 

Traschen, and R. Brandenberger, Phys. Rev. D51, 5438 (1995); L. Kofman, A. Linde, 

and A.A. Starobinsky, Phys. Rev. D56, 3258 (1997). 
[7] I.I. Tkachev, Phys. Lett. B376, 35 (1996); L. Kofman, A. Linde, and A.A. Starobinsky, 

Phys. Rev. Lett. 76, 1011 (1996); S. Khlebnikov, L. Kofman, A. Linde, and I. Tkachev, 

Phys. Rev. Lett. 81, 2012 (1998). 
[8] M.F. Parry and A.T. Sornborger, Phys. Rev. D60, 103504 (1999). 
[9] I.Tkachev, S. Khlebnikov, L. Kofman, and A. Linde, Phys. Lett. B440, 262 (1998). 
[10] S. Kasuya and M. Kawasaki, Phys. Rev. D58, 083516 (1998) 

[11] V. Zanchin, A. Maia, Jr., W. Craig, and R. Brandenberger, Phys. Rev. D60, 023505 
(1999). 

[12] L. Gammaitoni, et.al. "Stochastic Resonance", Rev. Mod. Phys. 70, 223 (1998) 
[13] S. Digal, R. Ray, S. Sengupta, and A.M. Srivastava, |hep-ph/9805227 . 
[14] A.M. Srivastava, Phys. Rev. D 47, 1324 (1993). 

[15] D. Boyanovsky, H.J. de Vega, R. Holman, and J.F.J. Salgado, Phys. Rev. D54, 7570 
(1996). 



7 



FIGURE CAPTIONS 



1) The central region is instantaneously heated to a temperature Tq — 0.72. A far 

separated vortex- ant ivort ex pair is produced due to flipping of $ in the heated region. 

2) $ at t = 908.0 for the resonant oscillation case for a portion of lattice, showing 
randomly oriented domains, with 4 vortex-antivortex pairs. Tq = 0.31 and Ta = 0.125 for 
this case. 

3) (a) Probability density of $ at t = 328.0 (b) Corresponding contour plot. Darker 
portions signify larger probability. 

4) Evolution of defect density n at large times. Straight lines show linear fits to the 
respective curves, n for a; = 1.19 and 1.31, are shown by the solid and the dotted curves 
respectively (with Tq = 0.5 and Ta = 0.125). Dashed curve shows n for Tq = 0.31, Ta — 
0.125 and u = 1.19. 

5) Variation of asymptotic average defect density with Tq. Dots show the values obtained 
from numerical simulations, while the curve shows Gaussian fit to these points. 
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FIG. 3. 
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